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We develop a quasi-chemical theory for the study of packing thermodynamics in dense liquids. 
The situation of hard-core interactions is addressed by considering the binding of solvent molecules 
to a precisely defined 'cavity' in order to assess the probability that the 'cavity' is entirely evacuated. 
The primitive quasi-chemical approximation corresponds to a extension of the Poisson distribution 
used as a default model in an information theory approach. This primitive quasi-chemical theory is in 
good qualitative agreement with the observations for the hard sphere fluid of occupancy distributions 
that are central to quasi-chemical theories but begins to be quantitatively erroneous for the equation 
of state in the dense liquid regime of pd 3 >0.6. How the quasi-chemical approach can be iterated to 
treat correlation effects is addressed. Consideration of neglected correlation effects leads to a simple 
model for the form of those contributions neglected by the primitive quasi-chemical approximation. 
These considerations, supported by simulation observations, identify a 'break away' phenomena 
that requires special thermodynamic consideration for the zero (0) occupancy case as distinct from 
the rest of the distribution. A empirical treatment leads to a one parameter model occupancy 
distribution that accurately fits the hard sphere equation of state and observed distributions. 



I. INTRODUCTION 

The quasi-chemical theory Jl], ^, || is a fresh attack on 
the molecular statistical thermodynamic theory of liq- 
uids. It is intended to be specifically appropriate in de- 
scribing liquids of genuinely chemical interest. But, in 
view of its generality, the quasi-chemical theory must be 
developed and tested for its description of the paradig- 
matic hard sphere fluid. In addition to the conceptual 
point, these developments are expected to be helpful in 
subsequent applications of the quasi-chemical theory to 
real solutions. 

The foundational virtues of the hard sphere fluid for 
the theory of liquids are widely recognizcdjPl and the 



be interrupted. Then we derive the new equation of state 
format, learn what we can by comparison of the primitive 
quasi-chemical approximation with Monte Carlo simula- 
tion results, study correlation contributions to propose 
an improved equation of state format, and finally exam- 
ine how this improved format works. 

Interestingly, though the some of these basic consider- 
ations are regarded as 'preliminary' Eqs. |^ or P, and the 



interest in this system continues to evolve [5|, |6|, |7], |8|. Re- 
cent developments of the theory of hydrophobic effects^ 
| 0, 0, g g, f§, in ad- 

dition the related quasi-chemical theory, have empha- 
sized again the significance of packing issues in a re- 
alistic molecular description of complex liquids. This 
paper studies the hard sphere fluid and develops de- 
fault models with utility in recent information theory 
approaches |, @ || Q || || || g| || || || . 

A compact derivation requires several preliminary re- 
sults, including brief specifications of the potential dis- 
tribution theorem, of the expression of chemical equi- 
librium, and of the quasi-chemical formulation. Addi- 
tionally, the notation here is not elsewhere standardized 
because these ideas are unconventional. The plan of the 
paper is thus to collect the necessary preliminary results 
in Appendix ^| so that the conceptual argument needn't 



formal identification of the equibrium ratio Eq. A4 , have 
been given before and have wide generality. 



II. A QUASI-CHEMICAL VIEW OF THE 
SOLVATION FREE ENERGY OF HARD CORE 
SOLUTES 

The preliminary results of Appendix |A| permit an 
attack on the solvation thermodynamics of hard core 
species built upon a simple device. Let's consider a so- 
lute A that doesn't interact with the solvent S molecules 
at all. We will consider formation of AS„ complexes and 
Fig. |l| depicts such a cluster. The interaction contribu- 
tion /3A/xa is zero and the quasi-chemical Eq. A5 ex- 
presses 



In 



= Inxo . 



(1) 



But the left side here is a test particle average for solute 
that rigidly excludes solvent molecules from the region 
defined by the indicator function b\j. If the region is 
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FIG. 1: An example of an AS7 cluster considered in the text. 
The nucleus (A) is visible in the center. Each of the ligands 
(S) overlaps the nucleus but no other ligand. 



TABLE I: 'Hit-or-Miss' Monte Carlo estimates |22j], as de- 
scribed in Appendix [b], of lnKn*- ' for unit diameter hard 
spheres and disks, respectively. Ki' ^^ (ir) and K2^°' = i§|- 

( 3 g The sample size was 24 G-configurations and the re- 
sults are believed to be accurate to the number of significant 
figures given. 
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spheres (3d) 


disks (2d) 


1 


1.43241 


1.14473 


2 


1.53915 


0.71321 


3 


0.56585 


-1.190 


4 


-1.4697 


-5.241 


5 


-4.684 


-13.77 


6 


-9.168 




7 


-15.46 





occupancy of the volume stenciled by 6aj- Thus the 

multiplicative factors of ps m Xn oc K n ^ps n are 
augmented by a self-consistent 'mean field' Ap8[. 

For reuse below, we summarize the technical results of 
this argument for hard core solutes, writing 



taken as defining a physically interesting molecule pair 
excluded volume, then the right side of Eq. |l| gives the 
negative of the excess chemical potential for the hard 
core solute defined by b\j ■ This is an example of the well 
known relation for hard core solutes In pp = — f3Apuc 
with 'HC denoting 'hard core'g |, [h| y, [jj, [l| vA 



IS, 13,171 IS, 19 



, |20J . This observation sheds light on the 
compensation of inner and outer sphere contributions to 
the quasi-chemical Eq. A5 but is not surprising. We then 
consider 'chemical' equilibria for binding of S molecules 
to the A molecule. Of course, there is no interaction be- 
tween the A molecule and the solvent molecules. The 
binding is just the occupancy by a solvent molecules of 
the 'cavity' defined by b^j- Combining these considera- 
tions gives 



HC 



hi 



£ 

m>l 



(2) 



The K m are well-de fined but typically computationally 
demanding; see Eq. A4. The evaluation of K m (°) will 
require few-body integrals over excluded volumes as is 
discussed in Appendix [b|. The primitive quasi-chemical 
approximation is 



HC 



In 



m>l 



(3) 



with A 
self-consistency 



a 'mean field' factor that achieves the 
condition £ n nK [ n ] p s n X n 

'ps n \ n . This amounts to an exten- 
sion of the Poisson distribution for use in an information 
theory procedure p8|. Here psK 1 °' ) = <n> is the expected 



; »o = r 



1 



E K mPS T 
m>l 



(4) 



This co mbin es Eq. y and the potential distribution theo- 
rem Eq. A2 for this problem. 



III. PRIMITIVE QUASI-CHEMICAL 
APPROXIMATION 



We can give a simple demonstration of the quantitative 
results of the primitive quasi-chemical theory by consid- 
ering the hard disk (2d) and hard sphere fluids (3d). Ta- 
ble I gives Monte Carlo estimates of the K,/ ) for those 
cases. The predicted distributions x n for two densities 
are in Figs. || and ^. Equation of state results f3Afi(p) for 
these systems predicted by this primitive quasi-chemical 
theory are shown in Figs. ^ and ^. The primitive quasi- 
chemical approximation is remarkably successfully in all 
qualitative respects, particularly in view of its simplic- 
ity. In particular, the predicted occupancy distributions 
such as shown in Fig. || are remarkably faithful to the 
data. Nevertheless, the equation of state predictions be- 
gin progressively to incur serious quantitative errors at 
liquid densities pd 3 > 0.6, (Fig. ||). 



IV. TEST OF THE EQUILIBRIUM RATIOS 

As a direct check on the primitive quasi-chemical 
mechanism, we can focus on testing ideal populations 



Eq. A12 as approximations to formally correct popula- 
tions Eq. A9. It is then natural to consider the ratios 



3 



In x 





FIG. 2: For the hard sphere fluid at psd 3 =0.277, comparison 
for n<5 of the Poisson distribution (solid curve) with primi- 
tive quasi-chemical distribution (dashed curve) implemented 
with the information theory constraint on the first moment 
nx n = 47rpgd 3 /3. The dots are the results of Monte Carlo 
simulation as discussed in Appendix |5J The primitive 
quasi-chemical default model depletes the probability of high- 
n and low-n constellations and enhances the probability near 
the mode. 

In x n 




FIG. 3: As in Fig. | but for p s d s =0.8. The error bars indi- 
cate the statistical uncertainty by showing the 67% confidence 
interval. 



FIG. 4: /3Ap(p) for the two dimensional hard disk fluid on the 
basis of the primitive quasi-chemical approximation (dashed 
lined). The Ree-Hoover 3,3 Pade approximant^] is the solid 
line and the dash-dot line is the first virial coefficient approx- 
imation. 




FIG. 5: /3Ap(p) for the three dimensional hard sphere fluid 
on the basis of the primitive quasi-chemical approximation 
(dashed line) . The solid line is the prediction of the Carnahan- 
Starling equation of state, taken as the accurate basis for 
comparison, and the dash-dot line is the first virial coefficient 
approximation. 



Xj/xQ — K n ps n . Consideration of these ratios corre- 
sponds to shifting the curves of Figs. || and || so that 
the initial point is at the common value (0,1). A specific 
example is shown in Fig. ^. Compared with this normal- 
ization, it is clear that the observed equilibrium ratios 
K„ are greater than the ideal ratios K„(°). 



V. CORRELATIONS 

A point of view here is that the geometric weighting 
with the A's of Eq. || establishes a mean field that adapts 
to the prescribed density. We now consider how to go 
beyond that mean field description. One idea is to ex- 
tract the features of the summand of Eq. ^ that would 



give purely geometric weighting and then to analyze what 
remains. To this end, we define £ = exp(/3A/i) and con- 
sult the formal identification of the equilibrium ratios 
Eqs. A4 . Thus we can rewrite Eq. ^ as 



C = i 



x °/ r ' 

m>l 



K, 



(0) 



C PS 



with 



X0/r 



-PAU T . 



(5) 



(6) 



The remarkable Eq. ^| is formally exact and hasn't been 
given before. The correlation factors x / m might, in prin- 
ciple, be investigated on the basis of simulation data and 
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FIG. 6: In 



(x n /x ) qca 
(x n /x ) sinl 



vs. n comparing for the hard sphere 



fluid the primitive quasi-chemical approximate populations 
with those observed by Monte Carlo simulation for pd 3 =0.8. 
This normalization focuses on the relative sizes of K n and 
K„' ' suggesting that K„ >K n ' ' even after the maxent 
reweighting. The variations are modest except for the change 
between n=0 and n=l. The error bars indicate the statistical 
uncertainty by showing the 67% confidence interval. In the 
middle of the distribution the statistical uncertainty is mostly 
from the uncertainty in the denominator factor xq. 



information theory analysis. That is likely to a special- 
ized nontrivial activity except of the lower density cases 
where the primitive quasi-chemical approximation is sat- 
isfactory. 



A. Iterating the Quasi-Chemical Analysis 

Alternatively, the quasi-chemical rules suggest natural 
theoretical appr oxim ation for the equilibrium ratios given 
formally by Eq. A4. Applying the rule Eq. ||, for n>0, 



= ((e-^)) Q K n ^C +1 



1+ E K m/nPS r 
rti>l 



n+l 



1 + E Km/- 
m>l 



(0) 



Ps' 



(7) 



The K m /„ can be understood by considering the chemical 
equilibrium 



AS,. Sr 



+ mS 



AS n S m , 



(8) 



i.e. the original AS„ cluster is the solute and it provides 
a nucleus for a constellation of m S particles, different in 
type for the S' species. How to address the calculation 
of the K m/n (°) is discussed in Appendix |c[ 

It is still helpful to focus on the populations even 
though more coefficients are involved now. To do this 
we consider 



-/3AU„ 



(9) 



m> 1 



and, to accomodate the additional factor of ( multiplying 
the terms m>l, rearrange so that 



C = i 



E <( e -' 3Al/ -)) ^ m ( ) P s m C m 

m> 1 

1- E ((e- 0AU -)) o Kj°)p s ^ 

m>l 



(10) 



This last equation is significant particularly because it 
suggests that a principal consequence of correlations can 
be a uniform reweighting of all coefficients m>l. Strik- 
ingly, that is exactly the suggestion of Fig. [| 

We can use this insight to push the argument further: 
the fact that the primitive quasi-chemical populations 
for m>l are correct relative to each other means that 
the quantities ((e _/3AC/m )) n are nearly exponentially de- 
pendent on m. For, in the first place, when the density 
is high, almost all the population is in the center of the 
distribution, and the Lagrange multipliers are negligibly 
affect by the relative reweighting of the m=0 term. Then 
the alteration of the original geometric weighting is liter- 
ally irrelevant. In the second place, when the density is 
sufficiently low, these correlation factors are nearly unity 
anyway. So we can accurately write 



c 



1 + A( PS ) KjVp s m \ m 



(11) 



m>l 



xo 'breaks away' from the rest of the distribution and 
requires individual consideration when the density is high 
enough that xo is sufficiently small due to correlation 
effects. Nevertheless 



c-i 

Co-1 



(12) 



where Co is the primitive quasi-chemical approximate 
value. Thus, when the primitive quasi-chemical approxi- 
mation is sufficiently accurate, the difficulty of evaluating 
the corrections should be much reduced. 

Though it would be interesting to calculate correlation 
corrections on the basis of Eq. [7] and Appendix a sim- 
pler, empirical approach suffices for our present purposes. 
This is because the discrepancies seen in Fig. [i] are sub- 
stantial but not problematic and, therefore, the required 
A(ps) is simple. In particular, the form 



(3 An = In 1 



(0) 



Ps 



m>l 



(13) 



conforms accurately to the Carnahan-Starling equation 
of state; see Fig. [?]. The literal coefficient in Eq. 13] was 



obtained by fitting on the basis of the Eq. 12 to minimize 
the discrepancy with the Carnahan-Starling equation of 
state. 

The occupancies predicted by this empirical model are 
depicted in Fig. |l 
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FIG. 7: Equation of state with the empirical correlation 
model, Eq. [13. (Lower panel) The solid line is the Carnahan- 
Starlingequation of state and the dashed line is the model 
of Eq. (Upper panel) Discrepancy: the empirical corre- 
lation model (Eq. |l3[ ) less the Carnahan-Starling value. The 
mean absolute discrepancy against Carnahan-Starling equa- 
tion of state is about 1% and the maximum discrepancy is 
less than 3%, nearly as good conformance to the Carnahan- 
Starling model as that model to simulation data. When the 
final empirical parameter was fitted using only psd 3 < 0.3, the 
mean absolute discrepancy hardly changed but the maximum 
discrepancy doubled. 




FIG. 8: Predicted occupancies with the empirical correlation 
model, Eq. for psd 3 =0.8. Compare to Fig. ^. Again the 
solid line is the Poisson distribution. 



VI. CONCLUDING DISCUSSION 

Our first goal was to work-out how the quasi-chemical 
theory, a fresh attack on the the statistical thermody- 
namic theory of liquids, applies to the paradigmatic hard 
sphere fluid. The second goal was to work-out theoret- 
ical approximation procedures that might assist in de- 
scribing dense liquids of non-spherical species. The new 
fundamental results here apply generally to 'hard core' 
molecular models. The primitive quasi-chemical approx- 
imation, the procedure for iterating the quasi-chemical 
analysis, and the recognition of the 'break away' phe- 
nomenon of Fig. [6] are likely to be helpful in understand- 
ing packing in dense molecular liquids, beyond the hard 
sphere fluid. For the hard sphere system specifically, we 
have obtained a simple occupancy model, Eq [l3[ that is 
likely to be helpful in a variety of other situations. 

One situation is the description of packing restrictions 
when the quasi-chemical theory is used to treat genuinely 
chemical interactions, for example in the study of hydra- 
tion of atomic ions in water 2j|. The issue of 'context 
hydrophobicity' associated with many molecular solutes, 
including molecular ions, in water can also be addressed 
on the basis of quasi-chemical calculations and the devel- 
opments here. 

Another situation of current interest is the theory 
of primitive hydrophobic effects that has recently been 
reborng g @ @g ^ (l§ g|, @ |l| g. 
An historical view has been that the initial issue of hy- 
drophobic effects was the hydration structures and ther- 
modynamics following from volume exclusion by non- 
polar molecules in liquid water. The balance of at- 
tractive forces that might produce drying phenomena 
was a secondary concern, except that 'drying' was al- 
ways present in the scaled particle models [£6|. With the 
convincing clarification of the first of these problems]^, 

| 0, § § H H HI © H H & the is- 

sue of drying phenomena has been now taken up more 
enthusiastically^, |27], [28]. In this context, we note that 
the striking success of the two-moment information mod- 
els and the Pratt-Chandler theory §, g, [k], [H], g [lj, [TJ, 



ia lft 17 18 m 20, 29 



is due, in part, to a fortuitous 
balance of a 'gaussian' approximation in the theory |18[ 
and a compensating disregard for drying possibilities [|30|; 
both of these compensating approximations are expected 
to be benign for small molecule solutes (3^]. One ingredi- 
ent in a better understanding of this situation is a careful 
solution to the case where drying phenomena are entirely 
absent. That ingredient is better in hand with the results 
above. 
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APPENDIX A: PRELIMINARY RESULTS 



1. Potential Distribution Theorem 



The potential distribution theorem]!], 0, 14, [5l] may 

be expressed as 



(Al) 



where p a is the density of molecules of type a (the 'so- 
lute' under consideration), z a — exp(/3^ CT ) is the absolute 
activity of that species, q CT is the single molecule partition 
function for that species, and V is the volume. The dou- 
ble brackets ((. . .)} indicate the average over the thermal 
motion of the solute and the solvent under the conditions 
of no interaction between them, and the averaged quan- 
tity is the Boltzmann factor of those interactions. The 
average indicated here is the ratio of the activity of an 
isolated solute, p a V/q ai divided by the absolute activity, 
z a , of the actual solute. Thus 



In 



(A2) 



This is a formal result to the extent that evaluation of 
the quantities on the right side typically will involve non- 
trivial calculations on many-body systems. 



2. Chemical Equilibrium 

The traditional chemical thermodynamic consideration 
of a chemical transformation such as 



nAA + n B B ^ n c C + n D D 



(A3) 



with the formal result of Eq. A2 leads to the formal ex- 
pression 



K 



PG nc PB nD 
PA nA pB nB 

(«^ A£/O »o 



-0AUn 



'/Li ' 

V , 



^A4) 



(((e-^A)) ^)n A ( {{e -/JA%)) ^) 

This should be compared to the textbook result for ideal 
gas systems That comparison shows that the sin- 
gle molecule partition functions are multiplicatively aug- 
mented by the test particle averages But otherwise 
the structure of this important result is unchanged. The 
conclusion here is that the equilibrium ratios are well- 
defined objects though formal to the extent that non- 
trivial computational effort would be required to evaluate 
them on the basis of molecular information. 



3. Quasi-chemical Theory 

The quasi-chemical develop starts from consideration 
of a distinguished molecule in the solution and seeks to 
evaluate the chemical potential on the basis of events 
occurring within a defined 'inner sphere.' For a species of 
type A, that definition is codified by specifying a function 
&Aj that is equal to one (1) when solution molecule j is 
inside the defined region and zero (0) otherwise. Our 
starting point can be 



(3Ap A = In a: 



- In 



,-0AU A 



, (A5) 



where xq is the fraction of A solute species with zero (0) 
neighbors in the defined region. AUa is the interaction 
energy of the solvent with the solute A that is treated 
as a test particle . T he potential distribution theory per- 
spective on Eq. JA5|) is 



j 

P -/3AU A 



n(i-M»c 



(A6) 



The first, or chemical term, of Eq. A5 can be analyzed 
with chemical concepts associated with the reactions 



AS n 



nS ^ AS,, 



(A7) 



Here the indicated complexes are composed of n solvent 
(S) molecules within the defined region. Remember that 
the A molecule is a 'distinguished' solute molecule consid- 
ered at the lowest non-zero concentration [Q. The frac- 
tional amount of A species with a given solvation number 
n can be described by a chemical equilibrium ratio 



PAS n 



PAS„ =0 PS' 



(A8) 



The p a are the number densities and, in particular, p$ 
is the bulk number density of solvent molecules since the 
distinguished A molecule is infinitely dilute. This nota- 
tion permits the normalized re-expression 



K nPs n 



1+ £ K mPs r 

rn>l 



Since this yields 



(A9) 



(A10) 



the original Eq. (A5) can be re-expressed as 



In 



E 

m>l 



A" 



.PS' 



7 



In 



-I3AU 



(All) 



The virtue of these rearrangements is that the natural 
first approximation is 



(0) _ 



1+ £ KjVps™' 



(A12) 



The K n ^ are equilibrium ratios for the chemical reac- 
tion (A7) in an ideal gas. This formulation and the ap- 
proximation of Eq. ( A12 ) are closely related Q to the 
quasi-chemical (or cluster-variation) approximations of 
Guggenheim |33j, Bethe Jm), and Kikuchi [j35| . 

This approach should have greatest utility where the 
chemical balances of Eq. A7 are dominated by inner 
sphere chemistry that can be captured with computa- 
tions on clusters. Such chemical interactions are often 
much larger than the outer sphere contribution, the right- 
most term of Eq. A5 . 

But that outer sphere contribution remains and can't 
be neglected forever. An interesting example based 
on simulation of liquid water was discussed recently 
There the x was estimated from molecular dynamics re- 
sults and the remainder, the outer-sphere contributions 
to /3A/i, were positive, suggesting domination of those 
outer-sphere contributions by the packing constraints 
studied here. A principal goal of the present work is the 
development of a reasonable approach for describing the 
packing issues necessary for treating those outer sphere 
contributions. 

Reiss and Merryj36| analyzed population relations 
analogous to Eq. A9 but with activities appearing in the 
place of densities and with coefficients, here the equilib- 
rium ratios K„, appropriately different. The additional 
formal point here is the replacement of the activity by 
the density that permits the identification of the K n in 
Eq. |^, and then further permits consideration of the mean 
field treatment Eq. [| on the basis of an information the- 
ory constraint when K„' ' will be used. At this stage, 
the quasi-chemical approximation achieves a particularly 
primitive character and deviates from the goal of bound- 
ing these thermodnamic quantities that was pursued by 
Reiss and Merry [j36|. 



APPENDIX B: CALCULATION OF THE K„ (0) 
FOR HARD SPHERES AND HARD DISKS 



The K„(°) sought for reaction Eq. A7 are given by 

*n (0) = ■ (Bl) 

(qs/V) qAS a= „ 



(See Eq. A4.) For this problem, q$ 



V/A s 3 with A s 
a thermal deBroglie wavelength for S but these momen- 
tum integrals cancel perfectly and are irrelevant as usual. 



Therefore, 



J d 3 n...J d 3 r n j f[ e(i,j 



(B2) 



j>i=l 



The notation f. d 3 rk indicates the three-dimensional 
spatial integral over the volume of the A-ball, a sphere 
of radius 1. The indicated integrand is thus 3n dimen- 
sional. The integrand is zero (0) if \ri — rj| < 1 (overlap) 
for any (ij) and one (1) otherwise. Thus the integral can 
be estimated by sampling n-point uniform placments in 
the A-ball and scoring the fraction of such placements 
that are free from overlaps between the n unit diameter 
S-spheres. This approach fails for n larger than those pre- 
sented in Table Q. But larger clusters were not observed 
in our simulation of the fluid, so our approach should be 
regarded as satisfactory. 

The analogous two dimensional procedure was used for 
the hard disk results. 



APPENDIX C: CALCULATION OF THE K 



m/n 



(") 



In contrast to Appendix [B], here the ratio sought is 

(0) _ 9AS^S m 



./r, 



(qs/V) QAS' 



(CI) 



corresponding to the reaction Eq. ^. Again, the explicit 
factors of V, the momentum integrals, and the factor of 
n! all cancel perfectly so that 



m\K (0 ] = I / d 3 n 

m/n \ f 1 
MS' 



m 

'I'm i II e (hj) 
AS', \J> i=1 



d 3 ~ 



(C2) 

Here the notation J AS , d 3 rk indicates an integral over 

the excluded volume of an AS^ complex to an S ligand. 
The SS excluded volume, the integrand, is the same as 
before. But the structure of the AS^ complex fluctuates 
and the volumes obtained for specific structures are aver- 
aged over these fluctuations. The brackets (. . .) indicate 
the average over the structures of the isolated AS^ com- 
plex. This averaging is permitted and governed by the 
non-trivial denominator that appears in Eq. |C1[ 

Operationally, the calculation can be much as in Ap- 
pendix ^ except for (a) averaging utlizing a Metropolis 
Monte Carlo calculation for the n ligand spheres in the 
star AS^; and (b) random placements of the m additional 
points are into a sphere of radius two (2) since that would 
fully enclose any conformation of the cluster. 



APPENDIX D: CALCULATION OF x„ FOR THE 
HARD SPHERE FLUID 

The probability that there are n — 1 points in a sphere 
of radius r, x„_i(r), can be obtained from the distri- 
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bution, 4:Trpsr 2 T> n (r), of the distance r to the nth near- 
est neighbor of an arbitrary point. The probability that 
there are no more than n — 1 molecules in the void is 
equal to the probability that the nth nearest neighbor is 
at least r away from the void center 



n-l 

E 

m=0 



x m (r) = iirpt 



D n (y)y 2 dy 



(Dl) 



Isobaric-isothermal Monte Carlo can be used to calcu- 
late T> n {r). x n (r) for a range of r can be obtained 
from the distributions T> n (r). To increase the accu- 
racy of the estimated V n {r) for rarely observed r, small 
and large, a specific point in the simulation volume was 
chosen, and the sampling probability was reweightcd 



by 



47r ( o s ri 2 e- 47r ' ,sri3 / 3 + C 



where r, is the distance 



from the chosen point to the jth nearest center for each 
configuration and C is an empirically chosen, dimensional 



constant. This importance sampling 



4irp s r 2 V n (r) 



( [4 7 rp s r 1 2 e- 47r ' ,sri3 / 3 + c] 6(r - r„)J> 



<[4 7 rp s r'i 2 e- 47r ^n 3 /3 + C']) 

3 (D2) 

is based upon the idea that (r) — e~ 47rpsr / 3 is 
the function |37j appropriate for a random distribution 
of spheres. This idea attempts to make the observed 
distribution of the distance to the nearest particle more 
nearly uniform. The constant C was included to avoid an 
unbounded weighting function. The denominator of Eq. 
D2 is just a normalizing factor on the distribution. The 
denominator of Eq. D2 merely provides a normalizing 
factor. 

Isobaric-isotermal ensembles of 108 and 256 hard 
spheres were sufficient. The Carnahan-Starling equation 
was used to find the /3p needed for a hard sphere simu- 
lation at each specific density. 
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